Application of Bo Chen, et al. 
Attorney Docket No. 5180-02-1-1 
Priority: U.S. Provisional Application No. 60/398,937 and 
Provisional Application No. 60/407,277 
Filed: July 26, 2002 and August 30, 2002 



METHOD FOR SPECTROPHOTOMETRIC 
BLOOD OXYGENATION MONITORING 



1 HEREBY CERTIFY THAT THIS PAPER OR FEE IS 
BEING DEPOSITED WITH THE UNITED STATES 
POSTAL SERVICE "EXPRESS MAIL POST OFFICE TO 
ADDRESSEE" SERVICE UNDER 37 C.F.R. 1.10 ON THE 
DATE INDICATED ABOVE AND IS ADDRESSED TO 
MAIL STOP PATENT APPLICATION ; COMMISSIONER 
OF PATENTS; P.O. BOX 1450, ALEXANDRIA, VA 22313-1450 

Elizabeth A. Dudek 



"EXPRESS MAIL" MAILING LABEL 

NUMBER EV 332038601 US 

DATE OF July 24, 2003 




1 



This invention was made with Government support under Contract No. 
1 R43NS045488-01 awarded by the Department of Health & Human Services. The 
Government has certain rights in the invention. 

Method for Spectrophotometric Blood Oxygenation Monitoring 

This application claims the benefit of the filing date of United States 
Provisional Applications 60/398,937, filed 26 July 2002, and 60/407,277 filed 30 
August 2002. 

BACKGROUND OF THE INVENTION 

1. Technical Field. 

[0001] This invention relates to methods for non-invasively determining 
biological tissue oxygenation in general, and to non-invasive methods utilizing 
near-infrared spectroscopy (NIRS) techniques in particular. 

2. Background Information. 

[0002] The molecule that carries the oxygen in the blood is hemoglobin. 

Oxygenated hemoglobin is called oxyhemoglobin (Hb0 2 ) and deoxygenated 
hemoglobin is deoxyhemoglobin (Hb). Total hemoglobin is the summation of the 
two states of hemoglobin (Total Hb = Hb0 2 + Hb), and is proportional to relative 
blood volume changes, provided that the hematocrit or hemoglobin concentration 
of the blood is unchanged. The mammalian cardiovascular system consists of a 
blood pumping mechanism (the heart), a blood transportation system (blood 
vessels), and a blood oxygenation system (the lungs). Blood oxygenated by the 
lungs passes through the heart and is pumped into the arterial vascular system. 
Under normal conditions, oxygenated arterial blood consists predominately of 
Hb02. Large arterial blood vessels branch off into smaller branches called 
arterioles, which profuse throughout biological tissue. The arterioles branch off 
into capillaries, the smallest blood vessels. In the capillaries, oxygen carried by 
hemoglobin is transported to the cells in the tissue, resulting in the release of 
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oxygen molecules (Hb0 2 => Hb). Under normal conditions, only a fraction of the 
Hb02 molecules give up oxygen to the tissue, depending on the cellular metabolic 
need. The capillaries then combine together into venuoles, the beginning of the 
venous circulatory system. Venuoles then combine into larger blood vessels 
called veins. The veins further combine and return to the heart, and then venous 
blood is pumped to the lungs. In the lungs, deoxygenated hemoglobin Hb collects 
oxygen becoming Hb02 again and the circulatory process is repeated. 
[0003] Oxygen saturation is defined as: 



In the arterial circulatory system under normal conditions, there is a high 
proportion of Hb0 2 to Hb, resulting in an arterial oxygen saturation (defined as 
Sa0 2 %) of 95-100%. After delivery of oxygen to tissue via the capillaries, the 
proportion of Hb0 2 to Hb decreases. Therefore, the measured oxygen saturation 
of venous blood (defined as Sv0 2 %) is lower and may be about 70%. 
[0004] One spectrophotometric method, called pulse oximetry, determines 

arterial oxygen saturation (Sa0 2 ) of peripheral tissue (i.e. finger, ear, nose) by 
monitoring pulsatile optical attenuation changes of detected light induced by 
pulsatile arterial blood volume changes in the arteriolar vascular system. The 
method of pulse oximetry requires pulsatile blood volume changes in order to 
make a measurement. Since venous blood is not pulsatile, pulse oximetry cannot 
provide any information about venous blood. 

[0005] Near-infrared spectroscopy (NIRS) is an optical spectrophotometric 

method of continually monitoring tissue oxygenation that does not require pulsatile 
blood volume to calculate parameters of clinical value. The NIRS method is 
based on the principle that light in the near-infrared range (700 to 1,000 nm) can 
pass easily through skin, bone and other tissues where it encounters hemoglobin 
located mainly within micro-circulation passages; e.g., capillaries, arterioles, and 
venuoles. Hemoglobin exposed to light in the near infra-red range has specific 
absorption spectra that varies depending on its oxidation state; i.e., 
oxyhemoglobin (Hb0 2 ) and deoxyhemoglobin (Hb) each act as a distinct 
chromophore. By using light sources that transmit near-infrared light at specific 
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different wavelengths, and measuring changes in transmitted or reflected light 
attenuation, concentration changes of the oxyhemoglobin (Hb0 2 ) and 
deoxyhemoglobin (Hb) can be monitored. The ability to continually monitor 
cerebral oxygenation levels is particularly valuable for those patients subject to a 
condition in which oxygenation levels in the brain may be compromised, leading to 
brain damage or death. 

[0006] The apparatus used in NIRS analysis typically includes a plurality of 

light sources, one or more light detectors for detecting reflected or transmitted 
light, and a processor for processing signals that represent the light emanating 
from the light source and the light detected by the light detector. Light sources 
such as light emitting diodes (LEDs) or laser diodes that produce light emissions 
in the wavelength range of 700-1 OOOnm at an intensity below that which would 
damage the biological tissue being examined are typically used. A photodiode or 
other light source detector is used to detect light reflected from or passed through 
the tissue being examined. The processor takes the signals from the light 
sources and the light detector and analyzes those signals in terms of their 
intensity and wave properties. 

[0007] It is known that relative changes of the concentrations of HbC>2 and 

Hb can be evaluated using apparatus similar to that described above, including a 
processor programmed to utilize a variant of the Beer-Lambert Law, which 
accounts for optical attenuation in a highly scattering medium like biological 
tissue. The modified Beer-Lambert Law can be expressed as: 

\-~ log(/// c ) x = a x * C * d * B x + G (Eqn2) 

wherein "A^" represents the optical attenuation in tissue at a particular wavelength 
X (units: optical density or OD); "l 0 " represents the incident light intensity (units: 
W/cm 2 ); "I" represents the detected light intensity; "a x " represents the wavelength 
dependent absorption coefficient of the chromophore (units: OD * cm" 1 * fxM" 1 ); "C" 
represents the concentration of chromophore (units: |iM); "d" represents the light 
source to detector (optode) separation distance (units: cm); "B x " represents the 
wavelength dependent light scattering differential pathlength factor (unitless); and 
"G" represents light attenuation due to scattering within tissue (units: OD). The 
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product of "d*B x " represents the effective pathlength of photon traveling through 
the tissue. 

[0008] Absolute measurement of chromophore concentration (C) is very 

difficult because G is unknown or difficult to ascertain. However, over a 
reasonable measuring period of several hours to days, G can be considered to 
remain constant, thereby allowing for the measurement of relative changes of 
chromophore from a zero reference baseline. Thus, if time ti marks the start of an 
optical measurement (i.e., a base line) and time t2 is an arbitrary point in time after 
ti, a change in attenuation (AA ) between ti and t2 can be calculated, and 
variables G and l 0 will cancel out providing that they remain constant. 
[0009] The change in chromophore concentration (AC = C(t 2 ) - C(ti)) can 

be determined from the change in attenuation AA, for example using the following 
equation derived from the modified Beer-Lambert Law: 

AA A = - log(/, 2 //„ ) A = a x * AC * d * B x (Eqn 3) 

Presently known NIRS algorithms that are designed to calculate the relative 
change in concentration of more than one chromophore use a multivariate form of 
Equation 2 or 3. To distinguish between, and to compute relative concentration 
changes in, oxyhemoglobin (AHb0 2 ) and deoxyhemoglobin (AHb), a minimum of 
two different wavelengths are typically used. The concentration of the Hb02and 
Hb within the examined tissue is determined in (jjnoles per liter of tissue (jllM). 
[0010] The above-described NIRS approach to determining oxygenation 

levels is useful, but it is limited in that it only provides information regarding a 
change in the level of oxygenation within the tissue. It does not provide a means 
for determining the absolute value of oxygen saturation within the biological 
tissue. 

[0011] At present, information regarding the relative contributions of 

venous and arterial blood within tissue examined by NIRS is either arbitrarily 
chosen or is determined by invasive sampling of the blood as a process 
independent from the NIRS examination. For example, It has been estimated that 
NIRS examined brain tissue comprising about 60 to 80% blood venous and about 
20 to 40% arterial blood. Blood samples from catheters placed in venous 
drainage sites such as the internal jugular vein, jugular bulb, or sagittal sinus have 
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been used to evaluate NIRS measurements. Results from animal studies have 
shown that NIRS interrogated tissue consists of a mixed vascular bed with a 
venous-to-arterial ratio of about 2:1 as determined from multiple linear regression 
analysis of sagittal sinus oxygen saturation (SSSO2) and arterial oxygen saturation 
(SaC>2). An expression representing the mixed venous / arterial oxygen saturation 
(SmvC>2) in NIRS examined tissue is shown by the equation: 

Srnv0 2 = Kv * Sv0 2 + Ka * Sa0 2 (Eqn.4) 

where "SvCV represents venous oxygen saturation; "SaCV represents arterial 
oxygen saturation; and Kv and Ka are the weighted venous and arterial 
contributions respectively, with Kv + Ka = 1 . The parameters Kv and Ka may 
have constant values, or they may be a function of SVO2 and SaC>2. Determined 
oxygen saturation from the internal jugular vein (SijvC>2), jugular bulb (Sjb02), or 
sagittal sinus (SssCb) can be used to represent SVO2. Therefore, the value of 
each term in Equation 4 is empirically determined, typically by discretely sampling 
or continuously monitoring and subsequently evaluating patient arterial and 
venous blood from tissue that the NIRS sensor is examining, and using regression 
analysis to determine the relative contributions of venous and arterial blood 
independent of the NIRS examination. 

[0012] To non-invasively determine oxygen saturation within tissue at 
certain depth, it is necessary to limit the influence from the superficial tissues. 
For example, to determine brain oxygen saturation of adult human with NIRS 
technology, the contamination from extracranial tissue (scalp and skull) must be 
eliminated or limited. 

[0013] What is needed, therefore, is a method for non-invasively 

determining the level of oxygen saturation within biological tissue that can 
determine the absolute oxygen saturation value rather than a change in level; a 
method that provides calibration means to account for energy losses (i.e. light 
attenuation) due to light scattering within tissue, other background absorption 
losses from biological compounds, and other unknown losses including measuring 
apparatus variability; and a method that can non-invasively determine oxygen 
saturation within tissue at certain depth by limiting the influence from the 
superficial tissues. 
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DISCLOSURE OF THE INVENTION 

[0014] It is, therefore, an object of the present invention to provide a 

method for non-invasively determining the absolute oxygen saturation value within 
biological tissue. 

[0015] It is a further object of the present invention to provide a method 

that provides calibration means to account for energy losses due to scattering as 
well as other background absorption from biological compounds. 
[0016] It is a still further object of the present invention to provide a method 

that can non-invasively determine oxygen saturation within tissue at certain depth 
that limits the influence from the superficial tissues. 

[0017] According to the present invention, a method and apparatus for non- 
invasively determining the blood oxygen saturation level within a subject's tissue 
is provided that utilizes a near infrared spectrophotometric (NIRS) sensor capable 
of transmitting a light signal into the tissue of a subject and sensing the light signal 
once it has passed through the tissue via transmittance or reflectance. The 
method includes the steps of: (1) transmitting a light signal into the subject's 
tissue, wherein the transmitted light signal includes a first wavelength, a second 
wavelength, and a third wavelength; (2) sensing a first intensity and a second 
intensity of the light signal, along the first, second, and third wavelengths after the 
light signal travels through the subject at a first and second predetermined 
distance; (3) determining an attenuation of the light signal for each of the first, 
second, and third wavelengths using the sensed first intensity and sensed second 
intensity of the first, second, and third wavelengths; (4) determining a difference in 
attenuation of the light signal between the first wavelength and the second 
wavelength, and between the first wavelength and the third wavelength; and (5) 
determining the blood oxygen saturation level within the subject's tissue using the 
difference in attenuation between the first wavelength and the second wavelength, 
and the difference in attenuation between the first wavelength and the third 
wavelength. 

[0018] The present method makes it possible to account for energy losses 
(i.e. light attenuation) due to light scattering within tissue, other background 
absorption losses from biological compounds, and other unknown losses including 
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measuring apparatus variability. By determining differential attenuation as a 
function of wavelength, the energy losses due to scattering as well as other 
background absorption from biological compounds are cancelled out or minimized 
relative to the attenuation attributable to deoxyhemoglobin, and attenuation 
attributable to oxyhemoglobin. 

[0019] In order to account for the resulting minimized differential 
attenuation attributable to tissue light scattering characteristics, fixed light 
absorbing components, and measuring apparatus characteristics, each of the 
parameters must be measured or calibrated out. Since direct measurement is 
difficult, calibration to empirically determined data combined with data developed 
using the NIRS sensor is performed by using regression techniques. The 
empirically determined data is collected at or about the same time the data is 
developed with the NIRS sensor. Once the calibration parameters associated 
with attenuation attributable to tissue light scattering characteristics, fixed light 
absorbing components, and measuring apparatus characteristics have been 
determined, the NIRS sensor can be calibrated. 

[0020] The calibrated sensor can then be used to accurately and non- 
invasive^ determine the total oxygen saturation level in the original subject tissue 
or other subject tissue. In addition, if effective pathlength of photon traveling 
through the tissue is known, for example, the separation distance ("d") between 
the light source to the light detector is known or is determinable, and the value of 
"B^, which represents the wavelength dependent light scattering differential 
pathlength factor, then the total amount of concentrations of deoxyhemoglobin 
(Hb) and oxyhemoglobin (Hb02) within the examined tissue can be determined 
using the present method and apparatus. 

[0021] The calibrated sensor can be used subsequently to calibrate similar 
sensors without having to invasively produce a blood sample. Hence, the present 
method and apparatus enables a non-invasive determination of the blood oxygen 
saturation level within tissue. For example, an operator can create reference 
values by sensing a light signal or other reference medium using the calibrated 
sensor. The operator can then calibrate an uncalibrated sensor by sensing the 
same light signal or reference medium, and subsequently adjusting the 
uncalibrated sensor into agreement with the calibrated sensor. Hence, once a 
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reference sensor is created, other similar sensors can be calibrated without the 
need for invasive procedure. 

[0022] There are, therefore, several advantages provided by the present 
method and apparatus. Those advantages include: 1) a practical non-invasive 
method and apparatus for determining oxygen saturation within tissue that can be 
used to determine the total blood oxygen saturation within tissue as opposed to a 
change in blood oxygen saturation; 2) a calibration method that accounts for 
energy iosses (e.g., Sight attenuation) due to light scattering within tissue, other 
background absorption losses from biological compounds, and other unknown 
losses including measuring apparatus variability; and 3) a practical non-invasive 
method and apparatus for determining oxygen saturation within tissue that can 
distinguish between the contribution of oxygen saturation attributable to venous 
blood and that saturation attributable to arterial blood; 4) a practical non-invasive 
method and apparatus for determining oxygen saturation within tissue at certain 
depth that limits the influence from the superficial tissues. 
[0023] In an alternative embodiment, aspects of the above-described 
methodology are combined with pulse oximetry techniques to provide a non- 
invasive method of distinguishing between blood oxygen saturation within tissue 
that is attributable to venous blood and that which is attributable to arterial blood. 
Pulse oximetry is used to determine arterial oxygen saturation, and the arterial 
oxygen saturation is, in turn, used to determine the venous oxygen saturation. 
[0024] These and other objects, features, and advantages of the present 
invention method and apparatus will become apparent in light of the detailed 
description of the invention provided below and the accompanying drawings. The 
methodology and apparatus described below constitute a preferred embodiment 
of the underlying invention and do not, therefore, constitute all aspects of the 
invention that will or may become apparent by one of skill in the art after 
consideration of the invention disclosed overall herein. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0025] FIG.1 is a diagrammatic representation of a NIRS sensor. 

[0026] FIG.2 is a diagrammatic representation of a NIRS sensor placed on 

a subject's head. 



9 

[0027] FIG. 3 is a diagrammatic view of a NIRS sensor. 

[0028] FIG.4 is a block diagram of the present methodology for calibrating a 

NIRS sensor. 

[0029] FIG. 5 is a graph showing an exemplary plot of absorption coefficient 

vs. wavelength. 

DETAILED DESCRIPTION THE INVENTION 

[0030] The present method of and apparatus for non-invasively determining 

the blood oxygen saturation level within a subject's tissue is provided that utilizes 
a near infrared spectrophotometric (NIRS) sensor that includes a transducer 
capable of transmitting a light signal into the tissue of a subject and sensing the 
light signal once it has passed through the tissue via transmittance or reflectance. 
The present method and apparatus can be used with a variety of NIRS sensors. 
The present method is not limited to use with this preferred NIRS sensor, 
however. 

[0031] Referring to FIGS. 1-5, the preferred NIRS sensor includes a 

transducer portion 10 and processor portion 12. The transducer portion 10 
includes an assembly housing 14 and a connector housing 16. The assembly 
housing 14, which is a flexible structure that can be attached directly to a subject's 
body, includes one or more light sources 18 and light detectors 19, 20. A 
disposable adhesive envelope or pad is used for mounting the assembly housing 
14 easily and securely to the subject's skin. Light signals of known but different 
wavelengths from the light sources 18 emit through a prism assembly 22. The 
light sources 18 are preferably laser diodes that emit light at a narrow spectral 
bandwidth at predetermined wavelengths. In one embodiment, the laser diodes 
are mounted within the connector housing 16. The laser diodes are optically 
interfaced with a fiber optic light guide to the prism assembly 22 that is disposed 
within the assembly housing 14. In a second embodiment, the light sources 18 
are mounted within the assembly housing 14. A first connector cable 26 connects 
the assembly housing 14 to the connector housing 16 and a second connector 
cable 28 connects the connector housing 16 to the processor portion 12. The 
light detector 20 includes one or more photodiodes. The photodiodes are also 
operably connected to the processor portion 12 via the first and second connector 
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cables 26, 28. The processor portion 12 includes a processor for processing light 
intensity signals from the light sources 18 and the light detectors 19, 20. 
[0032] The processor utilizes an algorithm that characterizes a change in 
attenuation as a function of the difference in attenuation between different 
wavelengths. The present method advantageously accounts for but minimizes 
the effects of pathlength and parameter "E", which represent energy losses (i.e. 
light attenuation) due to light scattering within tissue (G), other background 
absorption losses from biological compounds (F), and other unknown losses 
including measuring apparatus variability (N). E = G + F + IM. 
[0033] Refer to Figure 1 , the absorption A bx detected from the deep light 
detector 20 comprises of attenuation and energy loss from both the deep and 
shallow tissue, while the absorption A xx detected from the shallow light detector 19 
comprises of attenuation and energy loss from shallow tissue only. Absorptions 
A b?L and A xx can be expressed in the form of Equation 5 and Equation 6 below 
which is a modified version of Equation 2 that accounts for energy losses due to 
"E"" 

A u = -log(/ 1 ,// (> ) ;l = a A * C„ * L„ + a A * C, * L x + £ A (Eqn.5) 

Ki = ~ tog{U*. )a = «* * C * * L * + E ** (Eqn.6) 

Substituting Equation 6, into Equation 5 yields A\, which represents attenuation 
and energy loss from deep tissue only: 

4 =\x - Aa =«a * C * * L » -3a) 

(Eqn.7) 



= -log 



(It) 
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Where L is the effective pathlength of the photon traveling through the deep tissue 
and A'i and A'2 are the absorptions of two different wavelengths. Let 
E x = E x — E xX , therefore: 

A l —A 2 = AA 12 (Eqn . 8) 

Substituting Equation 7 into Equation 8 for A'i and A' 2 , AA'^can be expressed as: 

^12 = a m * C b* L b + ^#12 (Eqn.9) 
and rewritten Equation 9 in expanded form: 

A4 2 = ((a rl - a r2 )[Hb\ + - a o2 )[Hb0 2 \)L b + {E\ - E' 2 ) 
= ( Aa rl2 *[Hb\*L b ) + ( Aa ol2 *[Hb0 2 \ * L b ) + AE n 

(Eqn.10) 

where: 

represents the attenuation attributable to 

Hb; 

( ^ a o\2 * [Hb0 2 ]^ * L b ) represents the attenuation attributable to Hb0 2 ; 
and 

represents energy losses (i.e. light attenuation) due to light 
scattering within tissue, other background absorption losses from biological 
compounds, and other unknown losses including measuring apparatus 
variability. 
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[0034] The multivariate form of Equation 10 is used to determine [Hb0 2 ]b 
and [Hb] b with three different wavelengths: 



AA 12 - AE n 
AA n - AE' l3 



Aa rl2 Aa oU 
Aa rU Aa oU 

(Eqn. 11) 



[Hb\ 
[Hb0 2 \ 



Rearranging and solving for [Hb0 2 ]b and [Hb] b , simplifying the Aa matrix into 

[Act* ]: 



AA, 



12 



AA, 



13 



[Aa-r 1 ^)" 1 - 



AE 



12 



AE 



[Aa'J" 1 ^) -1 = 



13. 
(Eqn. 12) 



[Hb\ 
[Hb0 2 \ 



Then combined matrices [AA'] [Aa']" 1 = [Ac] and [AE] [AaT 1 = PP C ] : 



vi - 



l Hb 



v HbO, 



Hb 



HbO, 



[Hb0 2 \\ 



(Eqn. 13) 



The parameters A H b and A H b02 represent the product of the matrices [AAJ and 
[Aa 7 ]" 1 and the parameters W H b and W H b02 represent the product of the matrices 
[AE'J and [Act'] 1 - To determine the level of cerebral blood oxygen saturation 
(Sn02>, Equation 13 is rearranged using the form of Equation 1 and is expressed 
as follows: 



Sn0 2 % = 



[\b0 2 ^Hb0 2 ) 



[\b0 2 ^Hb0 2 + A m ^ Hb } 



*100% 



(Eqn.14) 
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Note that the effective pathlength U cancels out in the manipulation from 
Equation 13 to Equation 14. 



Equation 4 and the empirically determined values for SvC>2 and SaC>2. The 
empirically determined values for SvC>2 and Sa02 are based on data developed 
by discrete sampling or continuous monitoring of the subject's blood performed at 
or about the same time as the sensing of the tissue with the sensor. The temporal 
and physical proximity of the NIRS sensing and the development of the empirical 
data helps assure accuracy. The initial values for Kv and Ka within Equation 4 
are clinically reasonable values for the circumstances at hand. The values for 
A H b02 and AHb are determined mathematically using the values for \ bk and \ xx for 
each wavelength sensed with the NIRS sensor (e.g., using Equation 5 and 6). 
The calibration parameters W H b and ^Hb02, which account for energy losses due 
to scattering as well as other background absorption from biological compounds, 
are then determined using Equation 14 and non-linear regression techniques by 
correlation to different weighted values of SVO2 and Sa02; i.e., different values of 
Ka and Kv. Statistically acceptable values of Kv and Ka and ^Hb and ^Hb02 are 
converged upon using the non-linear regression techniques. Experimental 
findings show that after proper selection of Ka and Kv, the calibration parameters 
^Hb and WHb02 are constant within a statistically acceptable margin of error for an 
individual NIRS sensor used to monitor brain oxygenation on different human 
subjects. In other words, once the sensor is calibrated it can be used on various 
human subjects and produce accurate information for each human subject. The 
same is true for animal subjects. 

[0036] In an alternative method of determining the absolute oxygen 
saturation value Equation 7 is rewritten: 



[0035] 



The value for SnC>2 is initially determined from SmvC>2 using 




= a x * C * L b = {a rX [Hb\ + a oX [Hb0 2 \)L b 



(Eqn.15) 
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For a two wavelength system, let "R" be a calibration index parameter: 

R= Aj-ig; = (a rl [Hb] b + a ol [Hb0 2 \)L b 
Az-E', (a r2 [Hb\ + a o2 \Hb0 2 \)L b 



[HbQ 2 1 SnQ 2 

a +a \ Hb °A cc r2+ a o2 ^°^~ 
cc r2 +a o2 °*\-Sn0 2 



Canceling out U and substituting: 



[HbQ 1 \ _ SnQ 2 SnQ [HbQ 2 \ 

[Hb\ 1-Sn0 2 from 2 [Hb0 2 \+[Hb\ 



the following expression for SnC>2 is obtained: 

SnQ = a rl - a r2 R 

2 (a rl -a al ) + (a o2 -a r2 )R ^ n17 > 



[0037] The value of and A 2 ' are determined by measuring l b and l x for 
each wavelength. The parameters E'i and E' 2 can be considered as empirically 
determined calibration coefficients derived from the "best-fit" combinations of the 
weighted ratios of venous and arterial blood-oxygen saturation of the brain. By 
using non-linear regression techniques, the values of E'i and E'2 are determined 
by correlating to different combinations of venous and arterial oxygen saturation 
weighted values to find the "best-fit" relationship of "R" as a function of Ai\ A 2 \ ES 
and E'2 (Equation 17) to a specific ratio of venous and arterial saturation weighted 
values. 

[0038] In the determination of the SnC>2 percentage, the effective photon 
pathlength U cancels out. If, however, the photon pathlength is known or 
estimated, then the determination of the total value of Hb and/or Hb0 2 is possible. 
For example, if a value for pathlength L b is input into Equation 13 along with the 
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calibration values W H b and *F H b02, then the total value of Hb and/or Hb0 2 can be 
calculated. According to Equation 2, pathlength L can be estimated from the 
product of "B*d". The light source to detector separation (optode) distance 
parameter "d" in the pathlength calculation is a measurable value and can be 
made constant by setting a fixed distance between light source to detector in the 
NIRS sensor design. Alternatively, the parameter K d" can be measured once the 
optodes are placed on the subject by use of calipers, ruler, or other distance 
measurement means. The pathiength differential factor "B" is more difficult to 
measure and requires more sophisticated equipment. From a large data set of 
measured neonatal and adult head differential pathlength factor values, an 
estimation of the value of "B" can be determined within a statistically acceptable 
margin of error. Substitution of these predetermined values of "B" into Equation 13 
results in the determination of the total values of Hb and Hb0 2 . 
[0039] An alternative method of determining total values of Hb and Hb0 2 
combines Equation 3 and Equation 13 together. The multivariate form of 
Equation 3 is shown below: 



-10g(/, 2 //,lW L Al' 

log(/ f2 // tl ) A2 /L A2 
■log(/ r2 // (1 ) A3 /L A3 _ 



a HbX3 a Hb0 2 X3 



AHb 
AHbO, 



(Eqn.18) 



At time t = ti, the values of AHb and AHb0 2 are zero. Applying Equation 13, and 
knowing the calibration values of W hb and W H b02 at a predetermined differential 
pathlength factor "B" and optode separation "d", the total absolute values of Hb 
and Hb0 2 are determined at time t = ti, which are represented by [Hb] t i and 
[Hb0 2 ] ti respectively. At time t=t 2 , the values of AHb and AHb0 2 are then 
determined using Equation 18. The total values of Hb and Hb0 2 are then 
determined at time t = t 2 using the following equations: 



[Hb] =AHb(t 2 ) + [Hb] h 

(Eqn.19) 
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[Hb0 2 ] h = AHb0 2 (r 2 ) + [Hb0 2 ] h (Eqn . 20 ) 

Equations 19 and 20 are valid only if all the shared parameters in Equations 13 
and 18 are exact. Reduced to practice, the advantage of combining Equations 13 
and 18 result in improved signal to noise ratio (SNR) in the calculation of the total 
values for Hb and Hb0 2 . Conversely, improved SNR in the calculation of Sn0 2 is 
also obtained from the following expression: 



Sn0 2 %= , Hb ° 2 *100% , c ^ 

2 (Hb0 2 +Hb) ( Ec « n - 21 > 



[0040] After the calibration parameters Wnt and W H b02 are determined using 
the above-described methodology for an individual NIRS sensor, this particular 
sensor is said to be calibrated. A calibrated NIRS sensor affords accurate 
measurement of total tissue oxygen saturation, Sn0 2 , by non-invasive means. 
The calibrated sensor can be used thereafter on any human patient, including 
adults and neonates. The same is true for animal subject if the sensor was 
calibrated on animals. Although the present method is described above in terms 
of sensing blood oxygenation within cerebral tissue, the present method and 
apparatus are not limited to cerebral applications and can be used to determine 
blood oxygenation within tissue found elsewhere within the subject's body. 
[0041] According to an additional aspect of the present invention, the 
above-described method can also be used to establish a calibrated "reference" 
sensor that can be used to calibrate similar sensors through the use of a phantom 
sample (also referred to as a "reference sample"). The phantom sample has 
optical characteristics that are similar to the tissue being examined by the NIRS 
sensor. The calibrated reference NIRS sensor is used to sense the phantom 
sample and produce reference values. Similar, but uncalibrated, NIRS sensors 
can thereafter be calibrated by sensing the same phantom sample and adjusting 
either the hardware of the uncalibrated sensor or the output of the uncalibrated 
sensor until the output of the uncalibrated sensor agrees with the reference values 
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produced by the calibrated reference sensor. Therefore, the calibration 
parameters W H b and W H b02 for the uncalibrated sensor would be determined from 
the phantom sample. This technique makes it unnecessary to calibrate each new 
sensor in the manner described above, and thereby provides a relatively quick 
and cost effective way to calibrate NIRS sensors. 

[0042] Besides Hb and Hb02, other biological constituents of interest (e.g., 
cytochrome aa3, etc.) could be determined using the multivariate forms of 
equations 2, 3, 6 or 7. For each additional constituent to be determined, an 
additional measuring wavelength will be needed. 

[0043] In an alternative embodiment, the above-described methodology 
can be combined with pulse oximetry techniques to provide an alternative non- 
invasive method of distinguishing between oxygen saturation attributable to 
venous blood and that attributable to arterial blood. As demonstrated by 
Equation 4, SmvC>2 is determined by the ratio of venous oxygen saturation Sv0 2 
and arterial oxygen saturation Sa0 2 . A calibrated NIRS sensor affords accurate 
measurement of total tissue oxygen saturation, Sn02, by using regression 
techniques by correlation to mixed venous oxygen saturation Smv0 2 . Therefore, 
the following expression will result: 

Sn0 2 = Smv0 2 =K V * Sv0 2 +Ka* Sa0 2 (Eqn . 22) 

Non-invasive pulse oximetry techniques can be used to determine the arterial 
oxygen saturation (Sa0 2 ) of peripheral tissue (i.e. finger, ear, nose) by monitoring 
pulsatile optical attenuation changes of detected light induced by pulsatile arterial 
blood volume changes in the arteriolar vascular system. Arterial blood oxygen 
saturation determined by pulse oximetry is clinically denoted as Sp0 2 . If NIRS 
monitoring and pulse oximetry monitoring are done simultaneously and Sp0 2 is 
set equal to Sa0 2 in Equation 23, then venous oxygen saturation can be 
determined from the following expression: 
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Sv0 2 = 



Sn0 2 -(Ka*Sp0 2 ) 
_______ 



(Eqn. 23) 



For the brain, venous oxygen saturation Sv0 2 would be determined from internal 
jugular vein (Sijv0 2 ), jugular bulb (Sjb0 2 ), or sagittal sinus (SSSO2) and the 
parameters Ka and Kv would be empirically determined during the calibration of 
the NIRS sensor. Under most physiological conditions, Sp0 2 is representative of 
brain arterial oxygen saturation Sa0 2 . Therefore, depending on which venous 
saturation parameter was used to calibrate the NIRS sensor, this clinically 
important parameter (i.e., Sijv02, Sjb0 2 , or Sss0 2 ) can be determined by 
Equation 24 by non-invasive means. 

[0044] Since many changes and variations of the disclosed embodiment of 

the invention may be made without departing from the inventive concept, it is not 
intended to limit the invention otherwise than as required by the appended claims. 
[0045] What is claimed is: 



